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ABSTRACT 

A complete analytical and numerical treatment of all magnetohydro dynamic 
waves and instabilities for radially stratified, magnetized accretion disks is pre- 
sented. The instabilities are a possible source of anomalous transport. While 
recovering results on known hydrodynamic and both weak- and strong-field mag- 
netohydrodynamic perturbations, the full magnetohydrodynamic spectra for a 
realistic accretion disk model demonstrates a much richer variety of instabilities 
accessible to the plasma than previously realized. We show that both weakly 
and strongly magnetized accretion disks are prone to strong non-axisymmetric 
instabilities. The ability to characterize all waves arising in accretion disks holds 
great promise for magnetohydrodynamic spectroscopic analysis. 
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1. Hydromagnetic stability for stationary equilibria 

Magnetohydrodynamic (MHD) spectroscopy (Goedbloed et al. 1994) entails the ability 
to calculate all MHD perturbations accessible to a particular magnetized plasma configu- 
ration (the forward spectral problem), and holds the promise to use the MHD spectrum to 
diagnose the internal plasma state (the backward spectral problem). In this letter, we predict 
all MHD waves and instabilities for magnetized accretion disks. 

To analyse the entire MHD spectrum of gravitationally and thermally stratified, rotat- 
ing, magnetized equilibrium configurations, we use the formalism of Frieman and Rotenberg 
(1960). The equation of motion for the Lagrangian displacement £ of a fluid element is 

p^ + 2pv-V^ + Vn-B-VQ-Q-VB + V-(p£)g-V-[p£(v-V)v-pvv- V£] = 0, (1) 
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where II = — £• Vp— 7pV-£ + B-Q measures the Eulerian perturbation of total pressure, Q = 
V x (£ x B) is the Eulerian perturbation of the magnetic field B, and g is the gravitational 
acceleration. Equation (1) describes all waves supported by a time-invariant equilibrium 
with density p, pressure p, flow field v, gravity g, and B; 7 is the ratio of specific heats. We 
specify our analysis to one-dimensional (ID) axisymmetric stationary equilibria satisfying 



where the prime denotes differentiation with respect to the radial coordinate r from a 
(r, 9, z) cylindrical system. Hence, we concentrate on wave motions about an equilibrium 
where centrifugal forces, gravity, pressure gradients and magnetic forces balance. The one- 
dimensionality limits the generality by representing the gravitational influence of a central 
body of mass M* by a line source g(r) = —G c M*/r 2 e r , but this at most introduces a vertical 
wavelength cut-off on the validity of our results. 

The equilibrium symmetry allows for a decoupled analysis of individual normal mode 
solutions where £(r, 9, z, t) = (^ r (r),^g(r),^ z (r))expi(m9 + kz — u!t). For the radially varying 
magnetized equilibria, it is instructive to use a field line projection which introduces the basic 
variables x = r£ r , fj — ir (B/B x e r ) • £ and £ = irR/B ■ It is then a matter of algebra to 
turn the Frieman-Rotenberg system (1) into 



In this system, the radial derivative operator is denoted by D r , while the parallel gradient 
operator becomes the algebraic factor F = — iB ■ V = mBg/r + kB z . Furthermore, G = 
mB z /r — kB 9 and V g = vj/r — G c M*/r 2 . The latter measures the deviation from a Keplerian 
disk. The terms linearly proportional to the Doppler shifted frequency u = u — mvg jr — kv z 
represent the Coriolis effect. For static cylindrical equilibria without external gravitational 
field, this representation of the eigenvalue problem was introduced in Goedbloed (1975). 

From Equation (3), one can immediately obtain the Alfven and slow continuous parts of 
the ideal MHD spectrum by looking at extreme localization of the perturbations on single flux 
surfaces (D r — > 00). The first component can then be integrated at once, and, when inserted 
in the second and third, yields At) = and S( — 0, respectively. The Alfven continuum 
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corresponds to singular solutions fj oc S(r — ta), where A(ta) = (puo 2 — F 2 )(rA) = 0. The slow 
modes are given by ( oc 8(r — r s ), where S(r s ) = {pCb 2 {jp + B 2 ) — jpF 2 ) (r$) = 0. Due to 
the radial variation of the equilibrium and the Doppler shift, four continuous ranges of real 
eigenfrequencies are found. With the cluster points Q = ±oo for the fast subspectrum, these 
Doppler shifted forward and backward Alfven and slow continua determine the three-fold 
structure of the MHD spectrum. These continua are not influenced by gravity or rotation. 

An interesting limit of this system has obtained a lot of attention in the accretion disk 
literature. Considering a weakly magnetized accretion disk, it is possible to show that the 
combination of differential rotation with a weak magnetic field introduces a linear 'weak-field 
shearing' MHD instability (Balbus and Hawley 1991). As recently discussed in the review 
by Balbus and Hawley (1998), one can derive a sixth order dispersion relation governing local 
linear disturbances with purely vertical wave numbers k = ke z in disks where the restricted 
hydrodynamic equilibrium relation V g = holds. This dispersion relation for axisymmetric 
m = modes is formally found from Equation (3) by setting D r = 0, neglecting all curvature 
terms where Bg appears explicitly, and by using V g — 0. It can then be written as 
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Here, k 2 = 2vo(rvg)' /r 2 = G c M*jr z is the epicyclic frequency and the squared sound speed 
is c 2 = jp/p. This equation suggests marginal stability for F 2 /p = 'iGcM^/r 2 . 

An equivalent formulation of the MHD eigenvalue problem for equilibria satisfying (2) 
reduces the Frieman-Rotenberg system to a 2 x 2 system in terms of the variables (x, n): 
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In this equation, the following terms appear 
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This system (5) reduces to the normal mode picture by Appert et al. (1974) for static 
cylindrical equilibria without gravity. Bondeson et al. (1987) and Hameiri (1976) derived 
this system for cylindrical equilibria with flow, without gravity. An analysis for gravitating, 
flowing, magnetized equilibria in planar geometry is found in van der Hoist et al. (1999). 

The set of equations (5) again demonstrates that the continuous part of the MHD 
spectrum is found from A = and S = 0. The ranges in frequency where D = are not part 
of the continuous spectrum: for details, see Goedbloed (1998). Either of the formalisms given 
by the set (3), system (5), or the equivalent second order Sturm-Liouville type equation can 
be used to analyse the MHD spectrum of magnetized accretion disks. In particular, assuming 
a radial variation exp(ik r r) we obtain a local dispersion relation given by 



Similar to the static case, or the more general one including flow (Bondeson et al. 1987), one 
can prove the proportionality C 2 + DE oc AS, so that the continua AS = can be factored 
out. This leaves a sixth order polynomial in u> which governs all discrete local modes. Again, 
as a formal limit obtained by nullifying terms where m, k r , or Bg appears explicitly, and 
assuming a Keplerian equilibrium V g = 0, this polynomial is the dispersion relation (4). 
The more general relation (9) leads immediately to the following marginal stability criterion 
governing axisymmetric perturbations when Bg = 



The multiplicative factor F 2 = k 2 B 2 z highlights the essential magnetic character of this 
criterion. For non-vanishing Bg, even the m = case of equation (9) allows for overstable 
or damped MHD waves. The general relation (9) recovers various known results. Terquem 
and Papaloizou (1996) presented a stability analysis of accretion disks with both radial 
and vertical equilibrium variations, but a purely toroidal Bg magnetic field. In particular, 
a stability criterion for axisymmetric perturbations required a parameter defined in their 
equation (29) to be positive. If one considers only radial equilibrium variation, the same is 
found from our dispersion relation (9). Kim and Ostriker (2000) analyzed the stability of a 
cold (p = 0) radially stratified disk. Our relation (2) contains their equilibrium as a special 
case, and their results form a subset of our analysis. E.g., their equation (43) for 'poloidal 
buoyancy modes' follows from (9) by taking m = Bg = v z = and assuming k 2 B 2 /p 3> Co 2 . 
As a final note on the system (5), the pure hydrodynamical limit can be written as 



k 2 r ^ r + C 2 + DE = 0. 



(9) 
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The continuous part of the spectrum is now collapsed onto the flow continuum uj 2 = (Case 
1960). The local dispersion relation (9) with this flow continuum factored out reads 
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Several known results are immediately recovered from this relation, especially for m = 
modes where it has a simple solution. For a polytropic equilibrium with v z = 0, the limit 
where (k r ,k) — > 00 at a finite ratio k r /k yields uj 2 = n 2 k 2 /{k 2 + k 2 ), giving a dense range 
of discrete modes within —k < Co < k. The same limit shows that stability requires k 2 > 0, 
known as the Rayleigh (1916) criterion. The m 7^ global modes of equation (12) are 
of the Papaloizou and Pringle (1984) variety. High m non-axisymmetric instabilities can 
be inferred from (12), which simplifies in the incompressible limit (7 — > 00, p' — > 0) to a 
bi-squared relation. For weakly magnetized disks, the modes contained in the dispersion 
relation (12) will be complemented by Alfvenic and slow magnetosonic perturbations. 



2. Full spectral analysis of accretion disks 

With the general formalism for a complete normal mode analysis of magnetized accretion 
disks in place, we present in the remainder full MHD spectra for a specific thin disk model 
satisfying (2). We use an analytical model of a disk where the ratio (3 = 2p/B 2 is r- 
independent. Two other parameters are the constant helicity of the magnetic field a = 
—Bg/B z , and the disk aspect ratio e = H/r <C 1, where H is the disk scale height. If one 
neglects radiative pressure in the disk, e measures the ratio of the sound to the Keplerian 
speed e = c s /Qkr (Shakura and Sunyaev 1973). Because of the Newtonian gravitational 
field, we describe the physical quantities using radial power laws. To have a constant e, c s 
must scale as r -1 / 2 . The density profile is set to p oc r~ 3 / 2 . For the pressure and magnetic 
field we use the r~ l l 2 scaling of c s and the constant (3. The inner radius, the density and 
the Keplerian velocity at the inner radius make all quantities dimensionless. The resulting 
disk structure is: p = r~ 3 / 2 , p = e 2 r~ 5//2 , magnetic field B z = ey^r^ 5 / 4 / a/ (3(1 + a 2 ) and 
Bq = —aB z . The rotational velocity will then be v$ = v r~ l l 2 , where v a (a,(3,e) is found 
from (2). The v z component is set to zero. We calculate the collection of eigenfrequencies by 
a generalization of the LEDAFLOW program (Nijboer et al. 1997). LEDAFLOW computes 
the entire MHD spectrum of ID, gravitationally stratified, magnetized equilibria with flow 
by solving the linearised MHD equations. The radial discretization for the eigenfunctions 
employs quadratic and cubic Hermite finite elements, appropriate for both global and (near- 
singular local modes. We use rigid boundary conditions in the domain r e [l,r out = 10]. 
We input the parameters for the equilibrium (a, j3, e) and the normal mode numbers (m, k). 
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Since we consider thin disks, we set e — 0.1. The helicity of the B field is taken a — 1. 

Weakly magnetized accretion disk. We present two calculations for a weakly magnetized 
accretion disk. We set the parameter (3 = 2 10 3 . Similar to Balbus and Hawley (1998), 
we first take a purely vertical wave vector k = 2rr/He z . The vertical wavelength should 
not exceed the vertical size of the disk. In Fig. (1), we show the MHD spectrum. We can 
identify all modes in this plot: the fast-magnetosonic, epicyclic, and (overlapping) Alfven and 
slow-magnetosonic continuum modes are all on the real axis. A dense sequence of discrete 
modes appears on the imaginary axis: the magneto-rotational instability. This sequence is 
the unstable part of the slow magnetosonic subspectrum. Performing the same calculation 
for (3 — > oo (hydrodynamical case), only the epicyclic and sonic modes remain and no 
unstable mode is found, confirming its slow magnetosonic nature. We present in Fig. (2) the 
spectrum for non-axisymmetric perturbations with m — 10, k — 70. Note the asymmetry 
due to the Doppler shift k-v^0. The identification of the different modes is quite involved 
since the fast, slow, Alfven and epicyclic frequency ranges are all overlapping. Various 
branches of unstable modes can be seen, where only the rightmost branches correspond 
to their axisymmetric analogues already present in Figure (1). The structure suggests the 
presence of global modes accumulating towards both the Alfven and the slow continua. 
A full MHD spectroscopic analysis, which investigates the mode types from their specific 
polarization properties, is the subject of future work. 

Equipartition accretion disk ((3 = 1). We performed several calculations of MHD pertur- 
bations supported by disks where the magnetic pressure is of the same order as the thermal 
one. The results we have obtained are quite clear. If the perturbation is non-axisymmetric, 
unstable branches appear in the spectra. The growth rate of these instabilities grows for 
increasing m. For the disk equilibrium discussed above, no axisymmetric unstable modes 
are found. The spectral structure is shown in Fig. (3) for m = 10. Because all the terms in 
Eq. (5) are of the same order, a purely analytical approach is virtually impossible. 

3. Concluding remarks 

We applied a general formalism to an astrophysical system where most MHD instabil- 
ities are believed to occur: the accretion disk. Our analysis covers both unmagnetized and 
magnetized disks. Complete MHD spectra were presented for equilibria involving the gravity 
of a central object, velocities v$ and v z , and both thermal pressure and a magnetic field with 
azimuthal and vertical components. We deduced the relations (9), (10), (12), describing 
the threshold of several known instabilities for a differentially rotating fluid, including the 
Papaloizou-Pringle modes and the weak-field magneto-rotational one. The relations apply 
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to fully compressible, non-Keplerian magnetized disks. If one ignores vertical stratification, 
they generalize the analysis of Papaloizou and Szuszkiewicz (1992), who considered only 
m = modes when no azimuthal Bq field was present. We complemented our analytical 
approach with numerical calculations of MHD spectra of a realistic magnetized disk model. 
MHD spectra for weakly magnetized accretion disks contain the magneto-rotational insta- 
bility and several other branches of discrete unstable modes. We numerically found a variety 
of toroidal unstable modes which can affect equipartion accretion disks. Our results are con- 
sistent with Noguchi et al. (2000) who made an analysis of non-axisymmetric incompressible 
modes in accretion disks, and found Alfvenic instabilities for strong magnetic fields. We 
generalize their findings by considering all MHD modes. Instabilities at equipartition field 
strengths are of big interest since it is believed that magnetized disks launching MHD jets 
have a thermal pressure of the same order as the magnetic one (Ferreira and Pelletier 1995). 
Coppi and Coppi (1998, 2001) have pointed out the importance of anomalous transport 
produced by non-singular, non-axisymmetric bending modes for the regime of large mag- 
netic energy density. These modes are amongst the instabilities analysed here. A complete 
MHD spectroscopic study for these and other (also 2D) accretion disk equilibria will be pre- 
sented in future work (Goedbloed et al. 2002). Since the MHD spectrum for disks which are 
weakly magnetized or at equipartition field strengths shows this enormous variety of unstable 
modes, the description of MHD turbulence in magnetized accretion disks is far from com- 
plete. Follow-up nonlinear simulations should reveal the importance of the various modes for 
triggering and sustaining MHD turbulence. The analytical framework presented can be of 
interest for physical mechanisms involving disk instabilities, such as the 'Accretion-Ejection 
Instability' (Tagger and Pellat 1999) and the hydro dynamical Rossby vortices (Tagger 2001). 
The study of instabilities affecting radially stratified MHD jets (Kersale et al. 2000; Appl et 
al. 2000; Kim and Ostriker 2000) is embedded in the present formalism. 

This work was done under Euratom-FOM Association Agreement with financial support 
from NWO, Euratom, and the European Community's Human Potential Programme under 
contract HPRN-CT-2000-00153, PLATON, also acknowledged by F.C. 
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Fig. 1. — MHD spectrum for a weakly magnetized accretion disk for perturbations where k = 
ke z . As the Doppler shift vanishes (k.v = 0), the identification of the different subspectra is 
easy (A ± , S 11 * 1 for forward and backward Alfven and Slow modes, k for the epicyclic ones, and 
F* 1 for forward and backward fast modes.). The magneto- rotational instability is identified 
as a cluster spectrum of discrete modes associated with the slow continuum. 
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Fig. 2. — As Fig. (1), but for non-axisymmetric perturbations (m = 10). All the slow and 
Alfvenic sub-spectra overlap partially and several cluster branches of unstable modes are 
apparent. 
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Fig. 3. — MHD spectrum of toroidal perturbations in an equipartition accretion disk ((3 = 1). 
The calculation demonstrates a large collection of unstable modes. 



